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Abstract 

We study the spatial distribution of spin torque induced by a strong Rashba spin-orbit coupling 
(RSOC) in a ferromagnetic (FM) metal layer, using the Keldysh non-equilibrium Green's function 
method. In the presence of the s-d interaction between the non-equilibrium conduction electrons 
and the local magnetic moments, the RSOC effect induces a torque on the moments, which we 
term the Rashba spin torque. 

A correlation between the Rashba spin torque and the spatial spin current is presented in this 
work, clearly mapping the spatial distribution of Rashba spin torque in a nano-sized ferromag- 
netic device. When local magnetism is turned on, the out-of-plane {Sz) Spin Hall effect (SHE) 
is disrupted, but rather unexpectedly an in-plane (Sy) SHE is detected. We also study the effect 
of Rashba strength (a/j) and splitting exchange (A) on the non-equilibrium Rashba spin torque 
averaged over the device. Rashba spin torque allows an efficient transfer of spin momentum such 
that a typical switching field of 20 mT can be attained with a low current density of less than 
10'^A/cm2. 

PACS numbers: 72.25.-b,72.25.Ba,74.78.Na,75.75.-c 
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I. INTRODUCTION 



Ever since the theoretical prediction of the spin transfer torque (STT) [l, 2], there has 
been much research effort in utihzing the STT phenomenon to induce magnetization switch- 
ing and precession in ferromagnetic (FM) nanostructures without the need for an externally 
applied magnetic field. Devices which rely on the STT effect for magnetization switching 
offer the advantages of lower power consumption and reduced device dimension, which are 
crucial factors for nanoscale and high-density spintronic applications. The STT effect has 
been studied in conventional magnetic nanostructures such as spin valves jsl and magnetic 
tunneling junctions jj]. For STT to occur in these magnetic multilayers, one requires a pair 
of FM layers, i.e., a reference spin layer to generate a spin-polarized current for injection 
into the second free (switchable) layer. The two layers are magnetized in a noncollinear 
configuration so as to induce the transfer of the transverse spin momentum from the refer- 
ence to the free layer, which is mediated by conduction electrons flowing between the two 
layers. In the above process, the role of the spin-orbit coupling (SOC) effect is neglected. 
However, it is well-established that SOC can generate a nonequilibrium spin accumulation 
under the passage of current. Thus, it is conceivable that, in the presence of strong SOC 
effect, one can induce a STT without the need for an additional reference FM layer. This 
is corroborated by previous theoretical work which showed that the presence of Rashba 
spin-orbit coupling (RSOC) whose strength is denoted by a^j, and exchange interaction A 
between conduction electrons and local spins, can give rise to domain wall motion via spin 
momentum transfer The same spin transfer mechanism can also occur in a FM layer 
with a large and A values [g, 7|. The predicted RSOC-induced spin momentum transfer 
was experimentally demonstrated in a nanowire array a a. The above findings suggest 
that by utilizing Rashba-induced STT, one can achieve magnetization switching within a 
single FM layer, without an additional non-coUinear FM layer. Such single layer switching 
holds several potential advantages over conventional STT devices, such as a more symmet- 
ric current switching profile and the reduced influence of spin depolarization at the interfaces. 



A key element which determines the feasibility of Rashba STT is the presence of a strong 
Rashba SOC in the FM metal layer. Initial studies on the Rashba effect were focused on 
semiconductor (SC) materials 10l-ll4|. especially in two-dimensional electron gas (2DEG) 



hetero junction structures, which consist of two SC layers with different energy bandgaps. 
The conduction electrons in the 2DEG experience a strong RSOC effect due to the large 
potential gradient, as a result of the band-bending at the hetero junction interface. However, 
utilizing the Rashba-induced STT in SC materials is not an attractive proposition as SCs 
are intrinsically non-magnetic. Even if ferromagnetic behavior can be induced in them via 
doping (e.g. in dilute magnetic semiconductors or DMS), the resulting Curie temperature 
lies well below room temperature. Recent studies have shown, however, that a strong RSOC 
effect can also be induced in metallic nanostructures, both of the FM and non-FM types 
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18| . It is known that the Rashba SOC requires a structural inversion asymmetry (SIA), 



which gives rise to an internal electric field. In a metallic FM layer, the SIA can be enhanced 
by adjacent layers of heavy metals and oxides, which create the requisite band structure 
mismatch and large potential gradient at the interfaces 19|-|22|. By engineering the in- 



terfaces of the metallic FM layer, one can control the strength of the RSOC effect within 
the layer. The ability to enhance the RSOC coupling via interfacial effects has led to the 
experimental demonstration of the effect of Rashba-induced STT, as mentioned previously 
8|, |9[. However, to effectively harness this effect in future magnetic memory applications, it 
is essential to have an understanding of the microscopic spin transport in the presence of the 
RSOC effect, and the resulting non-equilibrium spatial distribution of the Rashba-induced 
STT. 



Thus, in this paper, we apply the Keldysh nonequilibrium Green's function (NEGF) 
technique to study the spin torque generated by the Rashba SOC on the local magnetization 
in a metallic FM layer. The NEGF method is suitable for the study of the Rashba STT, 
which is essentially driven by nonequilibrium spin accumulation generated by the passage 
of current in the presence of RSOC. In addition, the NEGF method can systematically 
incorporate the effects of the leads, and interactions (RSOC and exchange coupling) as self- 
energy terms. In Section II of the paper, we introduce the system Hamiltonian, consisting of 
the Rashba term Hso = oir{^ x p), where i; is a unit vector parallel to the internal electric 
field -E, which acts perpendicular to the FM layer, p is the electron wavevector and 
is the RSOC strength. The Hamiltonian also includes the s-d interaction characterized by 
the exchange energy A, which couples the nonequilibrium spin density due to RSOC effect 
to the local moments. Based on the second-quantized form of the Hamiltonian, we apply 
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the tight-binding NEGF formahsm, and calculate various microscopic transport quantities 
in the system, such as the local spin current and spin density, and the overall spin torque 
generated. In Section III, we numerically investigate (i) the spin torque efficiency as a 
function of the strengths of the RSOC effect {an), and the s-d exchange interaction (A), 

(ii) the relationship between the spin torque distribution and the local spin currents, and 

(iii) the in-plane spin Hall effect arising from RSOC. Finally, the summary of results and 
conclusion are presented in Section IV. 

II. THEORY AND MODEL 



FIG. 1. Schematic diagram of a ferromagnetic (FM) layer sandwiched between two dissimilar 
materials (oxides or heavy elements) to increase the vertical electric field and thus enhance the 
Rashba SOC effect. Current je flows in the in-plane x-direction. The magnetization of the FM 
layer M is oriented in the vertical z-direction. 

The structure under consideration is depicted in Fig. [H It consists of a metallic FM 
layer, sandwiched between two dissimilar materials (oxides or heavy elements) to enhance 
the RSOC interaction at the interfaces and within the FM layer. The local magnetization M 
is oriented along the vertical ^-direction. A charge current je is injected in the x-direction, 
which generates a field Hqq along the y = (zxje) direction. The Hamiltonian for the system 
can be expressed as: 



Z A 




H 



(1) 



P 



,2 




(2) 



(3) 
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where m is the free electron mass, and h is the reduced Planck's constant. Here, Hq denotes 
the kinetic energy of the conduction electrons in the FM layer, M is the magnetization 
direction, A is the exchange coupling between the free electron spin and the local moments, 
and S = {Sj} (where j = {x,y,z}) is the vector of Pauli spin matrices. Hgo denotes the 
Rashba interaction which couples the electron spin with its momentum, with p being the 
electron momentum, and the potential gradient inducing the RSOC effect being assumed 
to be in the direction z normal to the FM layer. The potential gradient may _arise from 

-3- 
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In 



a variety of sources such as impurities, host atoms, and structural confinement 
order to apply the many-body NEGF formalism, the above Hamiltonian has to be recast 
into the second quantized form: 



^0 = X] 4a(Cr±aa + Cr±ba) + ^raCl^Cra (4) 

r,(T ± r,(T 

Hso = —itsO ^ ^rai^i^x) aa'Cr±ba' T (Sy) ^^'Cr±a(T') ■ (5) 

where Cra{cl.„) is the fermionic annihilation(creation) operator of an electron with spin 
a =t, i at position r. Here, a is the lattice spacing representation on a square lattice in 
the tight-binding NEGF formulation , a = ae^ and b = aCy are the unit lattice vectors, to 
represents the hopping energy between lattice points, and is obtained by to = h/2'ma'^. The 
terms e^t = 4to + A/2 and = 4to — A/2 represent the on-site energy at the lattice site, 
and tso = o:R/2a is the SO coupling energy due to the Rashba interaction. 



In order to perform numerical analysis through the NEGF, the retarded (C) and lesser 
(G^) Green's functions are required. These are defined as 

G'r.yAt, t') = ^{{crAt), clAt')mt - t'), (6) 
G^.yAt,t')=^{clAt')crAt)) (7) 

After Fourier transformation, the expression for in energy space is given by 

G^{e) = [e-H-E^{e)]-\ (8) 
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In the above, S'' = is the retarded self-energy incurred by the lead a, where a = L{R) 

represents the left (right) lead. S'J^ can be determined by = Vag'^V^, where Va is the 
coupling matrix between the lead a and the FM layer, and is the retarded Green's function 
of the lead a and can be calculated numerically by the renormalization method 28||. 
can be calculated from the relation 

G'<(e) = G"'(e)S<(e)G"(e). (9) 

where = (G^^ . S< = S^, where = ifa^a is the lesser self-energy due to lead a, 
fa is the Fermi function in lead a, and Fq, = — 2ImEJJ^ is the linewidth function representing 
the coupling between the lead a and the central FM region. 



Various transport properties can be evaluated once the different Green's functions {G^ , 
G", and G*^) have been solved via Eqs. ([H]) and (jOj). The charge current through the system 
can be expressed in terms of the different Green's functions, as follows 

la =1 £^ de Tr {[S<(6)A(e)] - [T (e) aG< (e)]} , (10) 

where A{e) = ^[^^'(e) — G"(e)] is spectral function. Likewise, the current-driven local spin 
density in the central region is related to G^ as follows 

{Si)m = o ^ ^j i^i)cTCT' {clna-Cma') 



2 

o-cr' 



/CO 
deJ2iS^)..'G 
-oo __i 



mm,aa' \ ) ' 



oo 



where the subscript m refers to the site index. 



A / rfeTr[S,G<^(e)], (11) 



The spin torque exerted on the local magnetization can be defined as the difference 
between spin current going into and coiming out of the lattice point. We express the spin 
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torque as the divergence of spin current [29|: 



T = iXBjdVV-f\ (12) 

where V is the volume, /x_b is the Bohr magneton, and j^' is the spin current density between 
lattice points. The spin torque r can be also be defined as r = —'-/M x H, where 7 is 
the gyromagnetic ratio. We focus on the effective field H^g induced by Rashba SOC which 
acts on the local moments along the y = {z x j^,) direction. Thus, the effective field due to 
RSOC is 

^CfT = Hy = (13) 

where Mg is the saturation magnetization, and is obtained from Eq. (fT2|) . The torque 
efficiency is then given by ratio of Hcs/je- 



Under steady-state condition and in the absence of dissipative processes, the spin torque 
T, as defined according to Eq. ( IT2l) . is related to the divergence of the spin current. By 
considering the Heisenberg equation of motion, the local spin bond current between sites m 



and m' can be expressed in terms of 28|, |30|], i.e. 



(jSi \ ^ r-Si{kin)^ r.Si(SO)'. 
\Jmm'/ \Jmm' I \Jmm' /' 

02^^') = 'iEr^^ Tr[5,(G<,„(e) - G<„,(e))] 

0::::°^) = [e^ X (m' - m)]/^ J^^ 1^ Tr[(G<,^(e) + G<„,(e))], (14) 

where (m' — m) represents the unit vector between neighbouring sites on the x-y plane 
and Gj represents the unit vector of spin (Si). The above expression for the bond spin 
current comprises of two terms, i.e., the kinetic and SO coupling terms, arising from the 
corresponding terms in the Hamiltonian of Eqs. (|2]) and (|3]). By considering Eqs. ( !T2|) 
and ( IT4|) together, the spin torque is then given by the divergence of the spin bond current 
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(V ■ ji**), which in the discretized tight-binding model is approximated as:[28|, l3| 




(15) 



(16) 



(17) 



where Lso is the spin precession length (over which spin precesses by 1 radian), and can be 
expressed as L^o = '^^^ above constitutes to the spin torque expression of Eq. ([T2|) . 

III. RESULTS AND DISCUSSION 

Based on the tight-binding NEGF formulation presented in the above section, we per- 
formed numerical calculations of transport parameters such as the local spin density, bond 
spin current, and the effective field H^s in order to analyze the effect of RSOC induced non- 
equilibrium spatial spin torque on the FM layer structure. In our calculations, the following 
parameter values are assumed, unless otherwise stated: ur = 10^^^ eVm, m = 9.1 x 10^'^^ 
kg, a = 0.05 nm, Ep = 7.83 eV, M, = 1.09 x 10^ Am~\ A = 1.6 eV [31], and room 
temperature T = 300 K. The Fermi energy Ep and saturation magnetization Mg assume 
exemplary values corresponding to that of Co. 
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;j10° Acm"') 




;j10° Acm"') 




FIG. 2. The dependence of the effective current induced field (i^cff) due to the Rashba spin torque 
is plotted as a function of charge current density (je) for (a) varying Rashba strength an with a 
fixed exchange coupling A = 1.6 eV, and (b) varying exchange coupling A with a fixed a/j = lO"^'^ 
eVm. In (c), the spin torque efficiency (-ffeff/je) is plotted as a function of both A and a/j. In the 
calculations, we assume the dimension of the sample to be 50a x 50a, where a = 0.05 nm. 
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We first analyze the role of two key parameters aji and A in determining the strength of 



the effective field ifefr and the torque efficiency of the system. Figs. 2(a) and 2(b) show that, 



with a fixed and A respectively, H^s increases linearly with j^. This trend is consistent 
with the prediction that 



HeS = [Z X 



(18) 



derived from either gauge formulation \^ or from semiclassical (Boltzmann) transport equa- 
tion [7] in the strong coupling limit. Eq. ( ITS]) is a global expression of spin torque under 
linear response. In the gauge formulation, the factor P assumes a value of | in the adia- 
batic limit, while in the Boltzmann model, it refers to the spin polarization of current. We 
now consider the torque efficiency, which is given by the gradient of H^s with respect to 
je- As can be seen from Figs. 2(a) and 2(b), the torque efficiency is generally enhanced 



with increase in either an and A. However, in our non-equilibrium spatial treatment, it is 



clear from the plot in Fig. 2(c) that the torque efficiency does not vary linearly with an, 
unlike the prediction of Eq. ( JT8l) . The difference can be accounted for by noting that the 
global expression of Eq. ( [T8|) is derived in the limit of large coupling A, i.e., up to only 



the linear order in In our model, as can be seen from Fig. 2(c), the torque efficiency 



shows a slight oscillatory dependence superimposed upon a general increase with respect to 

an, especially at the region of < 10~^° eVm. However, at the region where > 10~^° 

eVm, its behavior is similar to the prediction derived from the Boltzmann semiclassical 

tl 

model for arbitrary coupling strength [7[. From the effective field Hes, one can estimate the 
critical current density required for magnetization switching. In Figs. 2(a) and 2(b), we 
consider RSOC strengths ranging from 10~^^ to 10"^'' eVm, which roughly corresponds to 
the practical values observed at the interfaces with heavy metal or oxide layers. Assuming 
an exemplary spin polarization of P = ~ 0.5, RSOC strength of oir = 10~^° eVm, and 
a switching field of Hs ~ 0.02 T applicable for Co nanowire structures [8|, we find that the 
critical current density for switching is approximately 10^ A/cm^ [see Fig. [2(^ . This IS 
significantly lower than the critical current density of the order of 10^ A/cm^ for the case of 



the conventional Slonczewski spin torque in spin valve structures 



32 



33| 
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5 10 15 20 25 30 35 40 45 50 



(a) 




(b) 

FIG. 3. The spatial distribution of the (a) Rashba effect spin torque Tx and its correlation with the 
local spin current by setting to 0.5 x 10"-'^'^ eVm, (b) and its correlation with (j^^/) 

by setting a/j to 1.5 x 10"^'^ eVm. The spin torque density is expressed in units of ^b/Lso- The 
sample has a lateral size of 50a x 50a. 
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Next, we examine the relationship between the Rashba- induced torque r and the spatial 



distribution of the spin currents. In Fig. 3(a), we plot the spin torque component based 
on the torque definition of Eq. fllSp . which relates it to the divergence of the local spin bond 
current jf!™/. For comparison, we plot the spatial distribution of the spin bond current j 



Sz 

mm' 



in Fig. 3(b) We observe a close correlation between the spatial distribution of and the 
fiow of the 2;-polarized spin current j^^/. The presence of RSOC causes a vortex-like flow 



of the bond spin current j^^, as shown in Fig. 3(b), Regions where j^^, is flowing in the 
+x {—x) direction corresponds to a large positve (negative) r^. Conversely, in regions where 
the positive and negative spin current fluxes meet and cancel each other, the spin torque 
becomes small. When the Rashba coupling strength an is increased, the magnitude of 



is generally larger since it scales with aji, as shown in Fig. 2(c) In addition, the vortices 
associated with the spin current become spatially smaller. This may be attributed to the 
increase in the rate of spin precession of the conduction electrons with an. The increased 
density of the vortices result in some cancelation of the bond spin currents near the center 
of the FM layer, so that more of the bond spin current flows at the boundaries, as shown in 
Fig. |3(b) 
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(c)A = 1.6 eV, = 5 X 10"" eVm, (d)A = 1.6 eV, UR^bx 10"" eVm, 

{sz)m- SHE disrupted {sy)m- SHE detecteded. 




5 10 15 20 25 30 35 40 45 50 



(e)A = 1.6 eV, ur = 1.5 x IQ-^'' eVm, 
{sy)m- SHE disrupted. 

FIG. 4. The spatial distribution of the spin density (a) {sz)m, (b) {sy)m, both with A = eV, 
aji = lx 10"i° eVm, (c) (s^)^, (d) {sy)m, both with A = 1.6 eV, or = 1.5 x 10"^° eVm. In (e) 
{sy)m is plotted with a larger aR = 1.5 x 10"^'' eVm, and A = 1.6 eV. The sample has a lateral 
size of 50a x 50a. 
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Finally, we analyze the spin density distribution and its dependence on the exchange 



strength A. Figs. 4(a) and 4(b) plot the spin density of {sz)m in the absence and presence 
of A, respectively. In the absence of exchange coupling (A = 0), the distribution profile of 
{sz)m clearly indicates a transverse separation of the 2;-spins, i.e. an out-of-plane spin Hall 
effect. This agrees with previous calculations based on the multimode scattering matrix 
method which predict a spin-Hall like separation of the out-of-plane spin component in 
the presence of Rashba effect [34]. However, the clear out-of-plane spin Hall separation 



disappears when a sizable exchange A is present, as shown in Fig. 4(c) It is found that the 
magnitude of {sz)m assumes a much larger value throughout the FM layer. This increase 
may be attributed to the alignment of the electron spin to the local moments oriented along 
the 2;-direction. We also analyze the in-plane spin density {sy)m distribution, as shown in 



Figs. 4(b) , 4(d) , and 4(e) There is no transverse separation of the in-plane spin density 
in the absence of A [Fig. 4(b)| . This is in line with theoretical prediction where the spin 
Hall effect induced by RSOC applies only to out-of-plane spins. However, in the presence 
of strong exchange coupling A, an "in-plane" spin Hall effect is present [Fig. 4(d)| . This 
in-plane spin Hall effect is destroyed in the presence of a strong Rashba strength, i.e. when 
aji is increased to 1.5 x 10^^° eVm [Fig. 4(e)| . This may be explained by noting that a 
large RSOC strength increases the rate of spin precession. Thus, the in-plane spin density 
{sy)m oscillates and changes signs along the direction of electron propagation (x-direction). 



as can be seen in Fig. 4(d) 



IV. CONCLUSION 

In summary, we have studied the non-equilibrium spatial intrinsic spin torque induced 
by Rashba spin orbit coupling in a ferromagnetic metal layer. Unlike the conventional 
Slonczewski spin torque, the Rashba induced torque is generated within a single layer, i.e. 
it does not require spin injection from an another ferromagnetic reference layer. We analyze 
the effect of two crucial parameters determining the strength of the Rashba spin torque: 
(i) the strength apt of the RSOC effect which is responsible for polarizing the injected 
charge current, and (ii) the exchange splitting A which couples the conduction electron to 
the local FM moments, thus allowing the transfer of spin momentum to the latter. The 
spin transport through the system is modeled via the tight-binding non-equilibrum Green's 
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function (NEGF) formalism. The NEGF theory systematically incorporates many-body 
effects including interactions with the leads as self-energy terms, and enables current and 
spin density to be evaluated spatially under nonequilibrium (bias-driven) conditions. Based 
on the NEGF theory, we numerically evaluate various transport parameters of the system, 
such as the effective field i^cff due to the spin torque, and the spatial distribution of the 
non-equilibrium spin current and spin accumulation. We found that H^s generally increases 
with both the RSOC strength ur and the exchange couphng A. However, the dependence 
of Hes on both parameters is not totally linear, unlike previous predictions based on gauge 
formulation or semiclassical Boltzmann which are global and only partially non-equilibrium 
(linear response), and in the strong coupling limits. For practical values of A and a^, the 
calculated critical current density corresponding to a typical switching field of 200 mT is 
calculated to be lower than 10^ A/cm^, comparable to that obtained via the conventional 
Slonczewski spin torque. For the structure under consideration where net current is in the 
x-direction and the local moments are aligned in the vertical 2;-direction, the net effective 
field (spin torque) is in the y (x)-direction. We plot the spatial profile of the x-component 
of the spin torque Tx, which bears a close correlation to that of the z-polarized bond spin 
current. It is also observed that the Rashba torque is concentrated near the boundaries 
of the FM layer. We also found that the combined presence of RSOC effect and exchange 
couphng A induces a Hall separation of in-plane spins, whereas the spin Hall effect for out- 
of- plane spins disappear with the introduction of A. Our calculations predict an effective 
field i^efif of the order of 1 Tesla for a current density of 10^ A/cm^, thus indicating the 
feasibility of utilizing the Rashba induced spin torque to achieve magnetization switching in 
spintronic applications. 
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